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The general methodology of binary switching requires tilting of potential landscape along the 
desired direction of switching. The tilt generates a torque along the direction of switching and the 
degree of tilt should be sufficient enough to beat thermal agitations with a tolerable error probability. 
However, we show here that such tilt is not necessary. Considering the full three-dimensional 
motion, we point out that the built-in dynamics can facilitate switching without requiring any 
asymmetry in potential landscape even in the presence of thermal noise. With experimentally 
feasible parameters, we theoretically demonstrate such intriguing possibility in electric field-induced 
magnetization switching of a magnetostrictive nanomagnet. 



Binary switching between two stable states of a 
bistable element is usually achieved by putting a bias field 
in the direction of desired state as conceived by Landauer 
and others |l|-l3| . The bias field lowers the potential pro- 
file along the desired state and thus creates an asymmetry 
in the energy profile. Another external agent inverts the 
potential profile making it monostable. It is necessary 
to make the monostable well deep enough to resist ther- 
mal fluctuations. Subsequently, the bias field is respon- 
sible for switching towards the desired direction when 
the potential landscape is turned back to bistable. This 
bias field whether it's a dc field or a field due to dipole 
coupling with an adjacent element has to overcome the 
thermal noise such that the switching takes place to the 
desired direction, however, with a permissible error prob- 
ability. Such widely accepted methodology is depicted in 
Fig. 

In this Letter, we show that the tilt or asymmetry in 
potential landscape is not necessary for switching, even in 
the presence of thermal noise. Our analysis considers the 
full three-dimensional space during a switching event and 
depicts the importance of considering the out-of-plane 
dynamics. We show that the built-in dynamics is suffi- 
cient to ensure error-resiliency. Such consequence cannot 
be unveiled without the consideration of out-of-plane dy- 
namics. All we need are a sufficiently high magnitude 
of the external agent inverting the potential landscape 
and a sufficiently fast ramp rate of the external agent. 
Basically, the internal dynamics of the system provides 
an equivalent asymmetry to cause error-resilient switch- 
ing. Such intriguing possibility is depicted in Fig. |l(b)| 
The potential energy landscape is never tilted to facili- 
tate switching in the desired direction. However, a torque 
persists due to out-of-plane motion, which ultimately di- 
rects switching in the desired direction. 

We theoretically demonstrate such feasibility in a mag- 
netostrictive single-domain nanomagnet, where stress is 
the external agent inverting the potential landscape of 
the nanomagnet. The magnetization of the nanomagnet 
switches successfully (> 99.99% probability) without any 
requirement of introducing barrier asymmetry even in the 
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FIG. 1. (a) The potential profile of a bistable switch is made 
monostable and tilted to switch from one state to another, (b) 
The potential profile remains symmetric but switching occurs 
due to built-in dynamics. 



presence of room-temperature thermal fluctuations. The 
stress can be generated on the magnetostrictive nano- 
magnet by applying a voltage on a piezoelectric layer elas- 
tically coupled to it. Such electric-field induced magne- 
tization switching mechanism has been studied recently 
because of their potential applications as extreme energy- 
efficient switches [4]. 

Switching dynamics of a single-domain magnetostric- 
tive particle in the presence of thermal fluctuations 
is modeled using the stochastic Landau-Lifshitz-Gilbert 
(LLG) equation [El-Q that describes the time-evolution 
of the magnetization under various torques. There are 
three torques to consider here: torque due to shape 
anisotropy of the nanomagnet, torque due to generated 
stress anisotropy, and the torque associated with random 
thermal fluctuations. Consider a magnetostrictive nano- 
magnet made of polycrystalline Terfenol-D that is shaped 
like an elliptical cylinder of major axis 100 nm, minor axis 
90 nm, and thickness 6 nm. These dimensions ensure that 
the nanomagnet has a single domain 



|. Fig. 2(a) shows 
the magnetization in three-dimensional space. We will 
call the z-axis the easy axis, the ?/-axis the in-plane hard 
axis, and the x-axis the out-of-plane hard axis based on 
the chosen dimensions of the nanomagnet. The two sta- 
ble magnetization states are along the ±z-axis. In order 
to switch the magnetization, uniaxial compressive stress 
is applied in the ^-direction. To analyze the magneti- 
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FIG. 2. (a) Magnetization in three-dimensional space, (b) 
The precessional motion due to applied stress tries to ro- 
tate the magnetization out-of-plane in anticlockwise direction 
while shape anisotropy torque due to small thickness of the 
nanomagnet tries to bring the magnetization to the magnet's 
plane (cj) = ±90°). Such counteraction happens only in the 
quadrants G (90°, 180°) and G (270°, 360°). 



zation dynamics, we will adopt the spherical coordinate 
system where the magnetization vector is in the radial 
direction, the polar angle is and the azimuthal angle is 
(j). Magnetization is on the magnet's plane if = ±90°. 

The total energy of the stressed nanomagnet is the sum 
of the shape anisotropy and stress anisotropy energy. The 
torque due to shape and stress anisotropy is derived from 
the gradient of potential profile and the random thermal 
torque is added separately in the stochastic LLG dynam- 
ics |7[. We will assume that the magnetization starts 
from 9 ~ 180° and the applied stress tries to switch it to 
9 ~ 0°. We will use the convention that magnetization's 
motion is in opposite direction to the torque exerted since 
the Lande ^-factor for electrons is negative. 

Dependence of shape-anisotropy energy on azimuthal 
angle <ft (rather than assuming <p — ±90°) generates 
an additional motion of magnetization, which is pro- 
portional to sin{2(j)) and vanishes when <j> = ±90°. 



As shown in the Fig. 2(b), the applied stress produces 
a torque that tries to rotate the magnetization anti- 
clockwise and forces it to reside out-of-plane. As the 
magnetization is deflected from the plane of the mag- 
net (0 = ±90°), the additional torque due to shape 
anisotropy as mentioned above [~sm(2^)] would try to 
bring the magnetization back to its plane. Because of 
such counteraction, out of the four quadrants for </>, i.e., 
(0°, 90°), (90°, 180°), (180°, 270°), and (270°, 360°), 
the magnetization would be stable in the second quad- 
rant or the fourth quadrant [i.e., (90°, 180°) or (270°, 
360°)]. Note that sin(20) is negative in these two quad- 
rants counteracting the precessional motion due to stress 
[see Fig. |2(b)] . We would call these two quadrants "good 
quadrants" and the other two (first and third) quadrants 
"bad quadrants", the reasoning behind which would be 
more prominent onwards. Consideration of the torque 
due to (^-dependence of shape anisotropy energy is in- 
strumental to the dynamics we present in this Letter. 
We will now describe the motion of magnetization in- 
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FIG. 3. Illustration of magnetization's motion in three- 
dimensional space under various torques generated due 
to shape and stress anisotropy alongwith considering the 
damping of magnetization (a is the phenomenological 
damping constant). Note that the dependence of shape 
anisotropy energy on <j> has generated two additional motions 
-\B shape ^{(j))\sin and -OL\B shape ^{(j))\ e<f,. The quadrant 
cj) G (90°, 180°) is chosen for illustration; choice of the other 
good quadrant cj) G (270°, 360°) is analogous. 



tuitively under various torques originating from shape 
and stress anisotropy as shown in the Fig. [3j We intend 
the motion of magnetization to be along the -eg direc- 
tion since we are switching magnetization from 9 ~ 180° 
towards 9 ~ 0°. The precessional motion of magneti- 
zation due to torque generated by the applied stress is 
in the -fe^ direction, but the damping of magnetiza- 
tion generates an additional motion, which is perpen- 
dicular to both the direction of magnetization (e r ) and 
-fe^, i.e., in -e# direction. These two motions are de- 
picted as 2B((j))cosO e<j> and —2aB((j))sin9cos9 e<9, respec- 
tively in Fig. |31 where a is the damping constant [Ej], 
B((j)) = B shape ((j))+B stress , BshapeW is the ^-dependent 
strength of shape anisotropy energy, and B s tress is the 
strength of the stress anisotropy energy Q . The quantity 
Bstress is negative and it must beat the shape anisotropy 
energy for switching to get started. Mathematically, note 
that both the quantities B((j>) and cos 9 are negative 
in the interval 180° > 9 > 90°. Hence, magnetization 
switches towards its desired direction due to the applied 
stress. However, this damped motion in -e# direction 
is weak because of the multiplicative factor, a, which is 
usually much less than one (e.g., a=0.1 for Terfenol-D). 

As magnetization rotates out-of-plane due to applied 
stress, and stays in the "good quadrants" for (j) [i.e., 
(90°, 180 °) or (270°, 360°)] as described earlier [see 
Fig. |2(b)| , it generates a motion of magnetization in 
the -ee (— x e r ) direction due to ^-dependence of 
shape anisotropy energy. Subsequently, a damped mo- 
tion is generated too in the -e^ (e r x — eg) direction. 
These two motions are depicted as — \B s hape,ci)( < i ) )\sin9 eg 
and -a\B s hape^((fy\ e</>, respectively in Fig. [3l where 
B shape, 0(0) ~ sin(2(j)). Note that in the good quadrants 
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FIG. 4. (a) Stress-cycle on the magnetostrictive nanomagnet, 

(b) Magnetization directions at different instants of time, and 

(c) Potential landscapes of the magnetostrictive nanomagnet 
in relaxed, compressively stressed, and expansively stressed 
conditions. 



for (/>, B s hape, (f)(0) is negative. Thus, keeping the mag- 
netization out-of-plane of the magnet in good quadrants 
is beneficial in switching the magnetization in its desired 
direction. In case the magnetization resides out-of-plane 
but in the bad quadrants, it would have resisted the mo- 
tion of magnetization in its desired direction of switching. 
A higher magnitude of stress keeps the magnetization 
more out-of-plane in good quadrants due to precessional 
motion, however, the damped motion —ot\B s hape,<f>(<t>) \ ^0 
tries to bring magnetization back towards the magnet's 
plane. As these two motions counteract each other, mag- 
netization keeps moving in the -eg direction and eventu- 
ally reaches at the x-y plane defined by = 90°. 

The stress-cycle alongwith the energy profiles and mag- 
netization directions at different instants of time is shown 
in the Fig. HI At time to, the magnetization direction is 
along the easy axis c± 180° and the potential landscape 
of the nanomagnet is unperturbed by stress. The poten- 
tial profile of the magnet is symmetric in both 0- and 
0-space with two degenerate minima at = 0°, 180° and 
a maximum at = 90° in #-space. The anisotropy in 
the barrier is due to shape anisotropy only, which is ~44 
kT at room-temperature using the magnet's dimensions 
and material parameters of Terfenol-D [7] . Note that the 
barrier height separating the two stable states (0 = 0° 
and 180°) is meant when the magnetization resides in- 
plane (i.e., <j) = ±90°) of the magnet. The barrier goes 



FIG. 5. Field and torque on the magnetization, M when it 
comes on the x-y plane (0 c± 90°). The field always tries to 
keep the magnetization in-plane (0 = ±90°) of the nanomag- 
net. Magnetization started from ~ 180°. (a) £ (0°,90°), 
(b) cj) £ (90°, 180°). Magnetization can traverse towards its 
destination c± 0° for the case (b), while it backtracks to- 
wards ~ 180° for the case (a). Note that the motion of 
magnetization is opposite to the direction of torque exerted 
on it. The quadrant <j> £ (180°, 270°) is analogous to the case 
(a), while the quadrant <j> £ (270°, 360°) is analogous to the 
case (b). 



higher when the magnetization is deflected from cj) = 90° 
as shown in the Fig. |4(c) [at time to]. The barrier is high- 
est when the magnetization points along the out-of-plane 
direction (cj) — 0° or 180°). Magnetization can start from 
any angle (\>i n iUai £ (0°,360°) in the presence of thermal 
noise [7[. 

As a compressive stress is ramped up on the nano- 
magnet between time instants to and £1, the potential 
landscape in #-space becomes monostable near <\> = ±90° 
provided a sufficient stress is applied. The potential bar- 
rier near <fi = 0° or 180° may not become monostable 
in #-space since barrier height is high therein, however, 
that is not necessary for switching. Since application 
of stress rotates the magnetization in ^-direction, it can 
eventually come near <fi = ±90° and starts switching from 
~ 180° towards = 90°. The minimum energy position 
between time instants t\ and £2 is at (0 — 90°, (j) = ±90°). 
From Fig. |4(c) , we can see that the potential profile at 
time instant t\ is still symmetric. 

Upon reaching at = 90°, if magnetization stays in the 
"good quadrants" for <j> [i.e., (90°, 180°) or (270°, 360°)], 
then the torque on the magnetization is in the correct 
direction so that it can traverse towards ~ 0°. (See 
Fig. [5j) This once again signifies the merit of termi- 
nology (good or bad) used for the four quadrants of 
(j). Note that at = 90° (i.e., cos0 = 0), the effect 
of stress on the magnetization rotation has diminished 
completely. The only two motions that are active at that 
point are -\B shape ^((j))\sin e e and -a\B shape ^((j))\ 
(see Fig. [3|). Since a <C 1, magnetization quickly gets 
out of the = 90° position and as the magnetization 
vector gets deflected from — 90° towards — 0°, the 
effect of stress comes into play. 

As the magnetization leaves from = 90° towards 
~ 0° and stress is ramped down at time £2 (see 
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FIG. 6. Percentage of successful switching events at room 
temperature (300 K) in a Terfenol-D/PZT multiferroic when 
subjected to stress between 10 MPa and 30 MPa. The critical 
stress at which switching becomes ^100% successful increases 
with ramp duration. However, at high ramp duration (e.g., 
120 ps), we may not achieve ^100% switching probability for 
any values of stress. 



Fig. H]), the torque due to stress tries to rotate the az- 
imuthal angle <fi of magnetization clockwise rather than 
anticlockwise (mathematically note that cosO is posi- 
tive for 90° > > 0° and B{(j)) is still negative when 
stress has not been brought down significantly, i.e., still 
\B s tress\ > \B shaV e(4>)\)- For a slow ramp-rate this ro- 
tation may be considerable and magnetization can stray 
into bad quadrants [(0°,90°) or (180°, 270°)]. Moreover, 
thermal fluctuations can aggravate the scenario. Switch- 
ing may impede and magnetization vector can backtrack 
towards where it started (0 ~ 180°) causing a switch- 
ing failure. In this way, switching failure may happen 
even after the magnetization has crossed the hard axis 
(0 = 90°) towards its destination (0 ~ 0°). This is why it 
does need a fast enough ramp rate during the ramp-down 
phase of stress. Note that when the potential landscape 
gets inverted from its unperturbed position (between the 
instants of time to and t\ in Fig. [4]), the effect of thermal 
fluctuations does not matter since it can only delay the 
magnetization to come at — 90°. 

Straying into a bad quadrant for azimuthal angle (j) 
[(0°,90°) or (180°, 270°)] during the ramp-down phase 
does not necessarily mean that the magnetization would 
fail to switch. Particularly if magnetization is close to its 
destination # ~ 0°, the shape anisotropy torque along- 
with the negative stress anisotropy torque (since stress is 
reversed) would have enough control to bring the magne- 
tization back to the other good quadrant for <j) and com- 
plete the switching. Thus, there may be ripples appear- 
ing in the magnetization dynamics at the end of switch- 
ing, which is because of the transition of azimuthal angle 
(j) between two good quadrants through one bad quad- 
rant [7]. 

Therefore, we require the following two criteria for suc- 
cessful switching without creating any asymmetry in the 
potential landscape even in the presence of thermal fluc- 
tuations: (1) a high enough stress that keeps the mag- 



netization more out-of-plane inside the good quadrants; 
and (2) a fast enough ramp rate that reduces the possibil- 
ity of backtracking of magnetization while it is crossing 
(or even after crossing) the hard axis (9 = 90°) towards 
its destination. A high stress and a fast ramp rate also 
increase the switching speed and surpass the detrimental 
effects of thermal fluctuations. 

Fig. [6] shows the switching probability for different val- 
ues of stress (10-30 MPa) and ramp durations (60 ps, 90 
ps, 120 ps) at room temperature (300 K). A moderately 
large number (10,000) of simulations were performed for 
each value of stress and ramp duration to generate these 
results. Initial angle distributions at 300 K for both 
and (j) are taken into account during simulations. The 
minimum stress needed to switch the magnetization at 
K is ~5 MPa, but the minimum stress needed to ensure 
switching at 300 K is ~14 MPa for 60 ps ramp duration 
and ~17 MPa for 90 ps ramp duration. For 120 ps ramp 
duration, ~100% success probability is unattainable for 
any values of stress since thermal fluctuations have higher 
latitude to divert the magnetization in wrong direction 
while stress is ramped down; at higher stress values ac- 
companied by a high ramp duration, there occurs higher 
out-of-plane excursion pushing the magnetization in bad 
quadrants, which further aggravates the error probabil- 
ity. 

In conclusion, we have shown that binary switching in 
a symmetric potential landscape is feasible, even in the 
presence of thermal noise. We have theoretically demon- 
strated such possibility in successful magnetization re- 
versal of a single magnetostrictive particle. The out- 
of-plane dynamics plays the crucial role to ensure such 
error-resilient switching. All we need are a sufficiently 
fast ramp rate and a sufficiently high magnitude of the 
stress applied on the magnetostrictive particle, which are 
experimentally feasible. 
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In this supplementary section, we provide further details of the analysis and some additional 
results. 

SI Magnetization dynamics of a magnetostrictive nanomagnet in 
the presence of thermal noise: Solution of the stochastic Landau- 
Lifshitz-Gilbert equation 

Consider an isolated nanomagnet in the shape of an elliptical cylinder whose elliptical cross section 
lies in the y-z plane with its major axis aligned along the z-direction and minor axis along the 
y-direction. (See Fig. [SlJ) The dimension of the major axis is a, that of the minor axis is 6, and 
the thickness is The volume of the nanomagnet is ft = (n/^abl. Let 6{t) be the angle subtended 
by the magnetization axis with the +z-axis at any instant of time t and <j)(t) be the angle between 
the +x-axis and the projection of the magnetization axis on the x-y plane. Thus, 9{t) is the polar 
angle and <j){t) is the azimuthal angle. Note that when <j) — ±90°, the magnetization vector lies in 
the plane of the magnet. Any deviation from <j) — ±90° corresponds to out-of-plane excursion. 

The total energy of the single-domain, magnetostrictive, polycrystalline nanomagnet, subjected 
to uniaxial stress along the easy axis (major axis of the ellipse) is the sum of the uniaxial shape 
anisotropy energy and the uniaxial stress anisotropy energy [S lj : 

E(t) =EsHA(t)+E S TA(t), (SI) 

where E$HA(t) is the uniaxial shape anisotropy energy and E$TA(t) is the uniaxial stress anisotropy 
energy at time t. The former is given by [SlJ 

E S HA{t) = (fi /2)M^N d (t) (S2) 

where M s is the saturation magnetization and N d (t) is the demagnetization factor expressed as [Slj 

N d (t) = N d _ zz cos 2 6{t) + N d _ yy sin 2 6(t) sin 2 cf)(t) + N d _ xx sin 2 6{t) cos 2 cf)(t) (S3) 

with N d _ zz , N d _ yyi and N d _ xx being the components of the demagnetization factor along the z-axis, 
y-axis, and x-axis, respectively. The parameters N d _ ZZ) N d _ yyi and N d - xx depend on the shape 
and dimensions of the nanomagnet and are determined from the prescription in Ref. [S2J. 

Note that uniaxial shape anisotropy will favor lining up the magnetization along the major axis 
(z-axis) by minimizing Esha{^)^ which is why we will call the major axis the "easy axis" and the 
minor axis (y-axis) the "in-plane hard axis" of the magnet. The x-axis will therefore be the "out-of- 
plane hard axis" of the magnet and it is "harder" than the in-plane one since the thickness is much 
smaller than the magnet's lateral dimensions (i.e., / << a, b). By constraining the nanomagnet from 
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Figure SI: An elliptical two-phase multiferroic nanomagnet stressed with an applied voltage. The 
voltage generates strain in the piezoelectric layer which is mostly transferred to the magnetostrictive 
layer. The resulting stress in the magnetostrictive layer favors aligning the magnetization vector in 
the plane defined by the in-plane hard axis (y-direction) and the out-of-plane hard axis (x-direction) 
rather than along the easy axis (z-direction) [S|3]. 



expanding or contracting in the ^/-direction, we can generate uniaxial stress along the z-axis (easy 
axis). An appropriate way to do this is to deposit the magnetostrictive layer on a piezoelectric 
layer to make a 2-phase multiferroic. By applying a potential across the piezoelectric layer, we 
generate a strain in it which is transferred to the magnetostrictive layer by elastic coupling. This 
generates uniaxial stress in the latter layer. It is the preferred technique for generating stress since 
it is electrical in nature and dissipates very little amount of energy [S3J. 
The stress anisotropy energy is given by [S lj 

Est A (t) = -(3/2)X 8 a(t)Q cos 2 6{t) (S4) 

where (3/2) A s is the magnetostriction coefficient of the nanomagnet and a(t) is the stress at an 
instant of time t. Note that a positive X s a(t) product will favor alignment of the magnetization along 
the major axis (z-axis), while a negative X s a(t) product will favor alignment along the minor axis 
(y-axis), because that will minimize Esta{^)- In our convention, a compressive stress is negative and 
tensile stress is positive. Therefore, in a material like Terfenol-D that has positive A 5 , a compressive 
stress will favor alignment along the minor axis, and tensile along the major axis. The situation 
will be opposite with nickel and cobalt that have negative X s . 

At any instant of time t, the total energy of the nanomagnet can be expressed as 

E(t) = E(0(t), <£(t), <r(t)) = B{t)sin 2 6(t) + C(t) (S5) 

where 

B(t) = B shape (t) + B stress (t) (S6a) 

B s ha P e(t) = B 8hape (c/>(t)) = yM s 2 fi [N d - xx cos 2 </>(t) + N d _ yy sin 2 c^{t) - N d _ zz ] (S6b) 
B s tress(t) = (3/2)X s a(t)ft (S6c) 

C(t) = ^M 2 QN d _ zz - (3/2)X 8 a(t)n. (S6d) 

Note that B stress (t) has the same sign as the X s a(t) product. It will be negative if we use stress to 
rotate the magnetization from the easy axis (z-direction) to the plane defined by the in-plane hard 
axis (y-direction) and the out-of-plane hard axis (x-direction) . 

In the macrospin approximation, the magnetization M(t) of the nanomagnet has a constant 
magnitude at any given temperature but a variable direction, so that we can represent it by the 
vector of unit norm n m (t) = M(t)/|M| = e r where e r is the unit vector in the radial direction in 
spherical coordinate system represented by (r, #,</>). This means that the magnetization vector is in 
the radial direction, and the polar and azimuthal angles will specify its orientation at 

any given time t. The unit vectors in the 9- and (^-directions are denoted by e# and e^, respectively. 
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The gradient of potential energy at any particular instant of time t is given by 

dE(t) 1 dE(t) 

VE(t) = VEWUMMt)) = ^e e+ -^y^e, (87) 



where 



and 



dE(t) 

— = 2B(t)sin9(t)cos6(t) (S8) 



= M*n(N d _ xx - N d _ yy )sin(2<P(t))sin z 9(t) 

= -B shape ^{t)sin 2 e{t). (S9) 



Here 



BshapeAt) = B shape A ( p(t)) = ^M*n(N d _ xx -N d _ yy )sin(2 ( j)(t)). (SfO) 

The torque acting on the magnetization per unit volume due to shape and stress anisotropy is 

T E (t) = -n m (t)xV£(fi(t)^(t),(7(t)) 

= -e r x [2B(t)sinO(t)cosO(t) e e - B shape ^(t)sin9(t) e^] 

= -2B(t)sin8(t)cos8(t) - B shape ^(t)sin6(t) e e , (Sll) 

whereas the torque due to thermal fluctuations is treated via a random magnetic field h(t) expressed 
as 

h(t) = Mt)e x + M*) 6 y + M*)e z (S12) 

where h x (t), h y (t), and h z (t) are the three components of h(t) in x-, y-, and z-direction, respectively 
in Cartesian coordinates. We will assume the following properties of the random field, h(t) [S4J. 

• The process h(t) is stationary. 

• The distribution of the quantities h x (t), h y (t), and is normal (Gaussian) with means 
equal to zero, i.e. (hi(t)) = where i = x,y, z. 

• The quantities hi(t) and hj(t f ) (where t 1 — t — ±A and i ^ j) are correlated only for time 
intervals A, which is much shorter than the time it takes for the magnetization vector to 
rotate appreciably. Furthermore 

(h i (t)h j (t)) = US ij 5(A) (i,j=x,y,z) (S13) 

where U = \^\(i+a^)M v [331 a * s ^ ne dimensionless phenomenological Gilbert damping 
constant, 7 = 2/j.B^o/h is the gyromagnetic ratio for electrons and is equal to 2.21 x 10 5 
(rad.m).(A.s) -1 , he is the Bohr magneton, and My — noM s Q. 

• The statistical properties of the quantities h x (t), h y (t), and h z (t) are isotropic. 
Accordingly, the random thermal field can be expressed as 



/ 2akT 

^ (t) = V | 7 |( 1 + ^)M v At G ^ (t) ( i = x >V> z ) ^ 

where 1/At is proportional to the attempt frequency of the random thermal field. Consequently, 
At should be the simulation time-step used to solve the coupled LLG equations numerically and 
G(o,i)(£) is a Gaussian distribution with zero mean and unit standard deviation [S6J. The simulation 
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time-step At should be selected small enough so that decreasing that step further does not make 

any significant difference in the results. Accordingly, the thermal torque can be written as 

T TH (t) = My n m (t) x h(t) = Pg(t) - P^t) e (S15) 

where 

p e (t) = My [h x (t) cos0(t) coscf)(t) + h y (t) cosO(t)sin</)(t) - h z (t) sin9(t)} (S16) 

P (t) = My [h y (t) coscf)(t) - h x (t) sin<f>(t)] . (S17) 

To derive the thermal torque, we have used the following identities. 

e x = sinO(t) cos<f){t) e r + cos9(t) cos<j){t) e# — sin(j){t) e^, (S18a) 

e y = sin6(t) sin(/)(t) e r + cos9(t) sin(/)(t) e# + cos<j){t) e^, (S18b) 

e z = cos9{t) e r — sin9{t) e#. (S18c) 



e r x e x = cos9(t) coscf)(t) + sin(f)(t) e#, (S19a) 
e r x e y = cos9(t) sin(j){t) — cos<j){t) e#, (S19b) 
e r x e z —sin9{t)e ( i ) . (S19c) 

The magnetization dynamics under the action of these two torques Te(£) and Tth(^) is de- 
scribed by the stochastic Landau-Lifshitz-Gilbert (LLG) equation as follows. 

dn m (t) / dn m (t)\ | 7 | 

-^-^ nm(t)X ^^J =-^[T E (t) + T TH (t)]. (S20) 

In the spherical coordinate system with constant magnitude of magnetization, the following 
relation holds: 

_ 5 _=_ e9 + „^(t)— e*. (S21) 

Accordingly, 

a ^n m (t) x dn ^*^ = -asin6{t) <f>'(t) e e + a9'{t) e (S22) 
where ()' denotes d{)/dt. This allows us to write 

dn ™^ - a ( n m (t) x = (9'(t) + asin0(t) <f/(t)) e e + (sm0(i) <f/(t) - a6'{t)) e . (S23) 

Equating the e# and components in both sides of Equation (jS20l) , we get 

0\t) + asind(t) <j)'{t) = -j^L {B shape ^(t)sin0(t) + P^t)) , (S24) 

My 

sinO(t) (j)'{t) - a0'(t) = -J^- (2B(t)sinO(t)cosO(t) - Pg(t)) . (S25) 

My 

Solving the above equations, we get the following coupled equations for the dynamics of 9{t) and 

(1 + a 2 ) ^ = ^- [B shape ^{t)sin9{t) - 2aB{t)sin6{t)cos6{t) + {aP e + P^)} . (S26) 

' (sinO + 0f)S27) 



(1+ ^<w> _ hi 



*B 3hapei< i,(t) + 2B(t)cOs9(t) - — (Pfl - ^) 
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S2 Initial fluctuation due to thermal torque 



We can see from Equation (IS2|) that because Nd- zz < ^d-yy <^ iV^_ ra , the shape anisotropy energy 
is minimum when 9 — 0° or 180°. Therefore, the magnetization of the unstressed magnet should 
be lined up along the easy axis (z-axis) in the absence of thermal perturbation. If that happens, 
then no amount of stress can budge the magnetization vector from this orientation since the torque 
due to stress vanishes when sinO — (see Equation (ISllj) ). Accordingly, 9 — 0° and 180° are called 
stagnation points since the magnetization vector stagnates at these locations. In this section, we 
show mathematically that thermal torque can overcome stagnation and deflect the magnetization 
vector from the easy axis {sinO — 0). 

When sinO — and no stress is applied on the nanomagnet, Equations (1S25|) and (1S24|) yield 

0'{t) = J^m) (S28) 

®0'(t) = j^Pe®- ( S29 ) 
Substituting for P$(t) and P^t) from Equations (jS16D and (jS17p . and using 9 = 180°, we get 

ah x (t)sin(f)(t) — ah y (t)cos(f)(t) = h x {t)cos(f){t) + h y {t)sin(j){t) (S30) 
which gives the expression for (f)(t): 

4>{t) = tan-' ( ak /*\ + } ! X ?l ) - (S31) 



ah x (t) - h y (t) 

If we substitute this value of <j)(t) in Equation fjS28D or in Equation fjS29D . we get 

t(ti = -|7l , "' (t> + "' (t> (S32) 

We can see clearly from the above equation that thermal torque can deflect the magnetization axis 
when it is exactly along the easy axis since the time rate of change of 0(t) is non-zero. Note that the 
initial deflection from the easy axis due to the thermal torque does not depend on the component 
of the random thermal field along the z-axis h z (t), which is a consequence of having ±z-axis as the 
easy axes of the nanomagnet. However, once the magnetization direction is even slightly deflected 
from the easy axis, all three components of the random thermal field along the x-, y-, and z-direction 
would come into play and affect the deflection. 

S3 Material parameters 

The material parameters that characterize the magnetostrictive layer are given in Table [STl [£(T1 SfTU] . 



Terfenol-D 


Major axis (a) 


100 nm 


Minor axis (b) 


90 nm 


Thickness (t) 


6 nm 


Young's modulus (Y) 


8xl0 10 Pa 


Magnetostrictive coefficient ((3/2)A s ) 


+90xl0- 5 


Saturation magnetization (M s ) 


8xl0 5 A/m 


Gilbert's damping constant (a) 


0.1 



Table SI: Material parameters for Terfenol-D as magnetostrictive layer. 



For the piezoelectric layer, we use lead-zirconate-titanate (PZT). The PZT layer is assumed to 
be four times thicker than the magnetostrictive layer so that any strain generated in it is transferred 
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almost completely to the magnet ostrictive layer. We will assume that the maximum strain that can 
be generated in the PZT layer is 500 ppm [SllJ, which would require a voltage of 111 mV because 
o?3i = 1.8e-10 m/V for PZT [S12J. The corresponding stress is the product of the generated strain 
(500 x 10 -6 ) and the Young's modulus of the magnetostrictive layer. So the maximum stress that 
can be generated on the Terfenol-D layer is 40 MPa. 



S4 Simulation results 

Some additional simulation results and corresponding discussions are given in Figures 



Distribution of e under thermal fluctuations 




165 170 
e (degree) 

(a) 



130 



Distribution of ^ under thermal fluctuations 




200 
e (degree) 

(b) 



400 



Figure S2: Distribution of polar angle ^initial an d azimuthal angle ^initial due to thermal fluctua- 
tions at room temperature (300 K). (a) Distribution of the polar angle 6 initial- The mean of the 
distribution is ^175°, while the most likely value is 180°. This is a nearly exponential distribution 
(Boltzmann-like) . (b) Distribution of the azimuthal angle (f) initial- These are two Gaussian distri- 
butions with peaks centered at 90° and 270° (or -90°), which means that the most likely location 
of the magnetization vector is in the plane of the nanomagnet. 
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Figure S3: Temporal evolution of polar angle 9{t) and azimuthal angle <j){t) for a fixed Qinmai — 175° 
and four different values of initial = {90°, 270°, 0°, 180°}. The applied stress is 15 MPa and the 
ramp duration is 60 ps. Thermal fluctuations have been ignored, (a) Temporal evolution when 
initial = 175°, ^initial = 90°. (b) Temporal evolution when 9 initia i = 175°, ^initial = 270°. (c) 
Temporal evolution when 9 initia i = 175°, ^initial = 0°. (d) Temporal evolution when 9 in iUai = 
175°, 6 

initial — 180°. Note that when 9 reaches 90° or even earlier, (j) always resides in the good 
quadrant [(90°, 180°) or (270°, 360°)], which makes the switching successful. For plots (a) and (b), 
the magnetization initially lies in the plane of the magnet ((^initial — ±90°) and precessional motion 
of magnetization due to applied stress is in the direction, which increases (f) with time. Thus, 
the magnetization starts out in the good quadrant. So both the motions (damped motion due to 
applied stress and the motion due to out-of-plane excursion depicted as —2aB((j))sin9cos9eo and 
— \B s } iape ^{(l))\sin9e0 ) respectively in the Fig. 3 of the main Letter) of magnetization are in the 
— eQ direction so that 9 decreases with time and the magnetization rotates in the correct direction 
towards 9 = 90°. The increasing out-of-plane excursion of the magnetization vector due to <j) 
increasing with time however is opposed by the damped motion due to out-of-plane excursion 
(depicted as — a\B s h ape ^{(j))\ in Fig. 3 of the main Letter), which tries to bring the magnetization 
back into the magnet's plane. These two effects quickly balance and <j) assumes a stable value in 
the good quadrant as seen in the plots (the flat regions of the 0-plots). When 9 reaches 90°, the 
torque due to stress and shape anisotropy vanishes. At this point, we start to reverse the stress 
and the damped motion due to stress and shape anisotropy eventually becomes again in the — e# 
direction. That continues to rotate the magnetization in the right direction towards = 0°, ending 
in successful switching. Slightly past 0.4 ns, continuing ^-rotation because of precessional motion 
due to stress and shape anisotropy pushes <j) into the bad quadrant, but eventually it escapes into 
the other good quadrant. This brief excursion into the bad quadrant causes the ripple in 9{t) versus 
t. For plots (c) and (d), the magnetization vector is initially lifted far out of the magnet's plane 
(^initial = 0°, 180°), where the huge out-of-plane shape anisotropy cannot be overcome by the stress 
anisotropy and B((f)) becomes positive, i.e., \B stress \ < \B s h ape ((j))\. So magnetization precesses in 
the clockwise direction (— e^) rather than in the anticlockwise direction (+e^). Thus <j) decreases 
with time, which immediately takes magnetization inside a good quadrant and eventually \B stress \ 
becomes greater than \B s h ape (<p)\. Then <j) assumes a stable value because of the damped motion due 
to out-of-plane excursion (depicted as —ot\B s h ape ^((j)) \ in Fig. 3 of the main Letter). Afterwards, 
switching occurs similarly as for the cases (a) and (b). 
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Distribution of delay [e=90 ] under thermal fluctuations 
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Distribution of § [e=90 ] under thermal fluctuations 
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Figure S4: Statistical distributions of different quantities when 15 MPa stress is applied with 60 
ps ramp duration on the magnet at room temperature (300 K). (a) Distribution of time taken 
for 9 to reach 90° starting from (6 initial A Hnitial) where the latter are picked from the distributions 
in the presence of thermal fluctuations (see Fig. IS2l) . The wide distribution is caused by: (1) 
the initial angle distributions in Fig. [S2] and (2) thermal fluctuations during the transition from 
some 9 = 0i n i t i a i to 90°. We point out here that we do need a sensing circuitry to detect when 
9 reaches around 90°, so that we can ramp down the stress thereafter. The sensing circuitry can 
be implemented with a spin-valve measurement of the magnetization. Such sensing circuitry is not 
uncommon in electronic circuits [£fT5]. Some tolerance is nonetheless required since the sensing 
circuitry cannot be perfect. Our simulation shows that the internal dynamics works correctly as 
long as the stress is ramped down when 9 is in the interval [85°, 110°], i.e. it does not have to 
be exactly 90°. If magnetization reaches at 9 — 90° (even past it) and stress is not withdrawn 
soon enough, then magnetization can end up at <j) — ±90° (potential energy minima), upon which 
the success probability would be 50% since thermal fluctuations can put it in either direction of 
the potential landscape. Thus, we take advantage of such non-equilibrium scenario for which a 
sufficiently high magnitude of stress and a sufficiently fast ramp rate are necessary, (b) Distribution 
of azimuthal angle (p when 9 reaches 90°. Note that <j) always resides in the good quadrant [(90°, 180°) 
or (270°, 360°)] and has a fairly narrow distribution. As explained in the main Letter, a high stress 
and fast ramp rate are required to ensure that (f) is in the good quadrants, which is conducive to 
successful switching. 
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Figure S5: Temporal evolution of the polar angle 9{t) and azimuthal angle <j)(t) when magnetization 
fails to switch and backtracks to the initial state. Simulations are carried out for room temperature 
(300 K). (a) The applied stress is 10 MPa and the ramp duration is 60 ps. When the polar angle 
9 reaches 90°, the azimuthal angle <j) has ventured into the bad quadrant (0°,90°). Thus, switching 
eventually fails, (b) The applied stress is 30 MPa and the ramp duration is 120 ps. When the polar 
angle 9 reaches 90°, the azimuthal angle <j) is greater than 90° and hence is in a good quadrant. 
However, after reaching 9 ~ 50°, the magnetization backtracks to the initial state and switching 
fails. This happens because of the long ramp duration and during this passage of time thermal 
fluctuations have ample opportunity to scuttle the switching by bringing <j) into the bad quadrant 
(0°,90°) 
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